Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 18 October 2012 



(MN style file v2.2) 



On the possibility of a warped disc origin of the inclined 
stellar discs at the Galactic Centre 



(N- 

A. Ulubay-Siddiki^*, H. Bartko^ and O. Gerhard^ 



o 
O 



< 

6 
in 



> 

o 

d 

T— I 

> 



^Istanbul University, Faculty of Science, Department of Physics, 34134 Vezneciler, Istanbul, Turkey 
^ Max-Planck-Institut fiir Extraterrestrische Physik,Giessenbachstrafle, D-85748 Garching, Germany 



18 October 2012 



1 INTRODUCTION 



The centre of our Galaxy hosts a supermassive black hole 



(iGenzel et alj 


20001: iGhez et al.l l2005l: iTripoe et all l2008l: 


Gillessen et al. 


20091). SgrA* is surrounded bv a cluster of 


old (TriDDe et all 20081: Schodel et al. |2009|). as well as a 


group of young stars ( 


Krabbe et al.ll 19951: Genzel et al.ll2003l: 


Levin & Beloborodov 


2003: Paumard et al. 20061: Lu et al. 


20061. 20091: Bartko et al.l 20091. 201(]f). 



Of the 136 young stars observed at distances 
pc to ~ 0.5 pc, 



Levin & Beloborodov' '20031; iPaumard et al.1 l2006l : iLu et al.) 
"2OO& . ,2009, : .Bartko ct al.. ,20lol ') which is observed to rotate 
clockwise (CW) on the plane of sky. Of the rest, 20 stars 
seem to populate an other disc highly inclined to the CW 
one (IGenzel et al. I l2003l : iPaumard et al. I I2OO6I : iBartko et al.l 
I2OO9I . I2OIOI )." and rotating counter clockwise (CCW) on the 



ABSTRACT 

The central parsec of our Galaxy hosts a population of young stars. At distances of 
r ^ 0.03 to 0.5 pc, most of these stars seem to form a system of mutually inclined discs 
of clockwise and counter clockwise rotating stars. We present a possible warped disc 
origin scenario for these stars assuming that an initially fiat accretion disc becomes 
warped due to a central radiation source via the Pringle instability, or due to a spinning 
black hole via the Bardeen-Petterson effect, before it cools, fragments, and forms stars. 
From simple arguments we show that this is plausible if the star formation efficiency is 
high, esF ^ 1, and the viscosity parameter a ~ 0.1. After fragmentation, we model the 
disc as a collection of concentric, circular rings tilted with respect to each other, and 
construct time evolution models of warped discs for mass ratios and other parameters 
relevant to the Galactic Centre environment, but for also more massive discs. We take 
into account the disc's self-gravity in the non-linear regime and the torques exerted 
by a slightly flattened surrounding star cluster. Our simulations show that a self- 
gravitating low-mass disc (Md/A/bh ^ 0.001) precesses with its integrity maintained 
in the life-time of the stars, but precesses essentially freely when the torques from 
a non-spherical cluster are included. An intermediate-mass disc (Afd/Afbh ~ 0.01) 
breaks into pieces which precess as independent discs in the self-gravity-only case, 
and become disrupted in the presence of the star cluster torques. Finally, for a high 
mass disc (Md/Mbh ~ 0.1) the evolution is dominated by self-gravity and the disc is 
broken but not dissolved. The time-scale after which the disc breaks into pieces scales 
almost linearly with Afd/A/bh for self-gravitating models. Typical values are longer 
than the age of the stars for Md/Afbh ^ 0.001, and are in the range ^ 8 x 10* — 10^ 
yr for Afd/A^bh 0.1 — 0.01 respectively. None of these discs explain the two Galactic 
Centre discs with their rotation properties. A comparison of the models with the 
better-defined clockwise rotating disc shows that the lowest mass model in a spherical 
star cluster matches the data best. 

sky (but also see lLu et al.l (|2006l )). Ages of these young stars 
are consistent with being a few Myr, suggesting that there 
has been a star formation episode in the Galactic Centre 
(GC) a few million years ago. 

In order for a molecular cloud near a SMBH to frag- 
ment into stars, its self-gravity should overcome the tidal 
field of the black hole. This requirement for star formation 
poses a constraint on the minimum cloud densities, which 
are orders of magnitude higher than the observed cloud 
densities near the GC. However theoretical estimates sug- 
gest that the fragmentation conditions are met naturally on 
the accretion discs which become self-gravitating beyond a 
few tenth of parsec ( Kolvkhalov fc Svunvaevlll980l : ICammij 



0.05 



59 popu l ate a disc (IGenzel et al. 20031: 
rodov' '2OO3I; IPaumard et al.ll2006l : ILu et aL 



I2OOII : iGoodmag 2qM )^ Therefore, several numerical simu- 
lations have been performed aiming at modeling the in-situ 
fragmentation of a nuclear/accretion disc for parameters rel- 
evant to the GC. The simulations were run either assum- 
ing a priori gravitationally unstable accretion disc already 
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in place (|Navakshin et all [200^ : [Alexander et aljliooj ). or 
trying to account also for the formation of the disc it- 
self through inf all of molecular cloud s into the vicinity of 
the black hole (iBonnell fc Ricd l2008l : iMapelli et al] l2008l : 
iHobbs fc Navakshinll2009f r 

Today it looks like a star forming disc at the GC can be 
simulated, albeit perhaps for a somewhat fine-tuned param- 
eter range. On the other hand, apart from the problem of 
youth, another issue still to be addressed is the distribution 
of the inclinations of the stars. It is reasonable to expect that 
a planar accretion disc leaves behind a planar distribution 
of stars when it frag ments, but the recent data published by 
iBartko et al.l (|2009l ) provide evidence for a warp in th e CW 
disc with an amplitude of about 60°. ILu et all (|2009l ) point 
out that even though the stars might have formed in-situ, 
their current orbital distribution suggests a more sophisti- 
cated origin t han a simple t hin a ccretion disc. Simulations 
performed bv lCuadra et al.l (|2008l ) are in line with this idea 
showing that once the stars form on a cold accretion disc 
it is not possible to perturb thes e stars to the high inclina - 
tions at which they are observed. iKocsis fc Tremaind (|201ll ) 
showed that vector resonant relaxation between the disc and 
the surrounding old star cluster might excite a warp in the 
disc if the stochastic torques from t he star clust e r dom inate 
over the disc's self-gravity. Recently, iHaas et al.l (|2010l ) con- 
sidered the evolution of the stellar disc using N-body sim- 
ulations. Their simulations successfully produced the rela- 
tive inclination between the stellar disc and the surrounding 
Circumnuclear Disc, however the origin of the disc was not 
addressed. 

Warped discs, although monitored only through maser 
emission from gas discs, exist on simila r scales in other 
nearby galactic n uclei such as NGC4258 |Herrnstein et al] 
Il996l). NGC1068 jGreenhiU fc GwinnI Il997l ). and Circinus 
I Greenhill et~aLl l2003h . An initially planar accretion disc 
could become warped when t o rqued by a spinning black hole 
llBardeen fc Pettersoiil Il975l: lArmitage fc NataraianI 1 19991 : 



iLodato fc PringldbOOTI : iMartirJ 200S ), when exposed to ra- 
diation from a central source "(|Pettersonlll977l : |Pringld[l99^ . 
I1997I ). or when subject to the torques induced by the stars in 
a ste lla r cusp around the black hole (|Bregman fc Alexander! 
I2OO9I ). iMilosavlievic fc Loebl (|2004l ') pointed out that the 
maser nuclei, and the Galactic Centre might represent dif- 
ferent epochs of a cycle during the lifetime of a typical spiral 
galaxy. 

In this paper, we investigate alternative scenarios for 
the formation of a star-forming, warped disc at the Galac- 
tic Centre. Accordingly, a flat accretion disc forms around 
SgrA*, extending out to the location of the young stars ob- 
served today. During a supposed period of active accretion, 
the disc is illuminated by the central source, or torqued 
by a spinning black hole, and becomes warped due to the 
Pringle instability or the Bardeen-Petterson effect, respec- 
tively. When the AGN activity subsides, the disc cools and 
forms stars. Afterwards, the stellar disc evolves in the grav- 
itational field of the black hole, its own self-gravity, and the 
surrounding old star cluster. We investigate under which 
conditions this scenario could work, and show that with a 
low mass such as inferred today the remnant warped stellar 
disc largely survives for the life-time of the observed young 
stars. 

In 32] we work out the conditions for which the disc 



would become warped due to the mechanisms mentioned 
above, and then fragment into stars. In !j3]we describe our 
model and the numerical scheme for studying the subsequent 
time evolution of the warped stellar disc. The results of our 
simulations are presented in i|4]where a comparison with the 
observations is also made. In JSJwe present a discussion, and 
in ^we summarize our results, and our conclusions. 



2 WARPING THE GALACTIC CENTRE DISC 

In this section, we discuss a plausible scenario of how a 
disc of young stars in the Galactic Centre might have ac- 
quired its warped shape. We start with the assumption 
that an accretion disc builds up, leading to an active phase 
of sub-Eddington accretion onto the Galactic Centre black 
hole. We then investigate two possible me c hanisms : radi- 
ation pressure instability (jPettersonl Il977l: iPringld Il996l. 
1997 ) and Bardeen-Petterson effect (jBardeen fc PettersonI 
for warping the accretion disc, which have both 



19751) 



been extensively discussed in the context of the maser 



Scheuer fc Feileill 19961: iPringk 


I1997I: lArmitage & NataraianI 


I999I: Lodato fc Pringld 2007: 


Martin 2008h. 



We further assume that after some time the accretion 
and energy production is reduced, which is followed by a pe- 
riod where the warped disc can cool and form stars. There- 
after the stellar disc is only subject to gravity, precessing 
under the influence of the gravitational torques from the 
disc itself, and the surrounding old star cluster. 

In H2.ll we constrain the surface density of the disc prior 
to fragmentation from the observed number density of young 
stars. In the following subsections 12.21 and 12.31 we consider 
in turn warping by the radiation pressure instability and the 
Bardeen-Petterson effect. In i ]2.4l we compare the radiation, 
viscous, and gravitational torques on the disc, and in i)2.5l 
we consider fragmentation and star formation. 



2.1 Surface Density of the Disc Prior to 
Fragmentation 

The warping mechanisms which we will briefly consider in 
sections (|2.2|) and (|2.3|l are generally studied within the 
framework of viscous, steady-state accretion discs. As such 
they make use of a number of parameters which from the 
observations of the stellar discs can not be tested or con- 
strained. For these, we will mostly refer to the canonical 
values for AGN discs when needed. Still, there is one pa- 
rameter, the surface density of the supposed gaseous disc 
which may be determined by the observations. To do so we 
make use of the stellar number counts on the GC discs. The 
total number of stars A'^* in a stellar population can be cal- 
culated by writing 



rM2 

m = ^{M)dM, 

J Ml 



(1) 



where ^(M) is the initial mass function, IMF, and Mi and 
M2 are the lowest and highest stellar masses assumed to 
exist in this population. The IMF describes how the mass is 
distributed in stars in a newly born population. For a given 
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Figure 1. Surface density of the supposed GC gaseous disc for 
star formation efficiencies of 0.01 (solid line), and 1 (dotted line). 



IMF, the total mass in stars, M. , is calculated from 



Mi{M)dM. 



(2) 



iBartko et all (|2010l ) d educe an IMF of i (M) = ^oA'/""-*^ 
for the GC discs. The lBartko et all (|2010l ) sample includes 
59 stars in the CW disc, and 20 stars in the CCW disc. 
Using equations ([T]) and and assuming a lower mass 
end of IMq, and an upper mass end of 12QMq, the cur- 
rent total stellar mass in the discs can be found to be 
M,|cw + Ms I CCW = Ms ~ 536OM0, where M^lcw, and 
Mslccw are the masses of the CW and CCW rotating discs 
respectively. The mass of the Galactic central black hole is 
~ 4 X 10® Mq, so the inferred disc mass corresponds to a 
mass ratio Md/Mbh = 0.00134. 

The total stellar disc mass inferred today is a fraction 
of the mass of the original gaseous disc since presumably 
not all the gas was converted to stars. Assuming a star for- 
mation efficiency esp the total mass of the seed disc can be 
calculated. Since the radial extent of the discs is observation- 
ally constrained, this mass can be converted into a surface 
density E of the proposed gas disc. The observations of the 
stellar dis cs suggest that the mass density decreases nearly 
as 1/r^ * jBartko et al-llioiol ). so we write for the mass Md 
of the gaseous disc 

Md = 5360^5360 espM©, (3) 

where /^536o = (Ms/536OM0) and for the surface density 
Ed(r) = So.if-i-'* we find 

So.i = 9.52 X loVsseoesF^ = 19.89^*5360 egp-^, (4) 
pc^ cm^ 

with f = r/O.lpc. In Fig. [1] we show for two star formation 
efficiencies the surface density of the disc obtained as de- 
scribed above. The solid line assumes esF = 0.01, and the 
dotted line assumes esF = 1- 



2.2 Radiation Driven Warping and the Galactic 
Centre Disc 

Radiation warp ing of accretion discs is studied in d etail by 
many authors (jPringld Il996l : iMalonev etal] 1 19961 : iPringld 



1 19971 : lOgilvie fc Dubusll200ll ). When an optically thick disc 
is exposed to central radiation, and it re-emits the absorbed 
incident photons parallel to the disc local normal, an inward 
directed force is experienced by each side of the disc. If the 
disc is slightly distorted, a net torque is induced which re- 
sults in a modification of the warp, and precession of the disc 
around the total angular momentum direction. Whether the 
disc will acquire a pronounced warp or not depends on the 
competition between the net torque, and the component of 
the stress in the r — z direction in cylindrical symmetry, 
present in warped discs. The latter forces the disc to set- 
tle onto a plane on an alignment time-scale. The alignment 
time scale at distance r from the black hole is associated 
with the r — z stress, and can be written as t,,^ = 2r^/i^2, 
where 1^2 is the vertical viscosity coefficient. The condition 
that the warp growth time scale, t r, be smaller than the 
viscous time scale, t^^ is written as (|Pringlelll997h 

127rEr^f7c 2r^ 

where L is the luminosity of the central source, c is the 
speed of light, and a is t he Shakura-Sunyaev parameter 
(jShakura fc SunvaevI 1197^ ). Equation ([5]) assumes an ac- 
cretion disc which is in steady state. For such discs, the 
luminosity of the disc L is related to the radiative effi- 
ciency e = L/Mc^, where M is the mass accretion given 
by A'l — 37rEi/i, and vi is the radial viscosity coefficient. 
In order to evaluate the warping criterion given in ((SJ one 
has to estimate the magnitude of the vertical viscosity 1/2. 
Previous studies of linear hydrodynamic warps and magne- 
tized shearing box simu lations of ac cretion discs showed that 
V2/V1 = 1/2q^ (IPapaloizou fc Prinaic 1983: Ogil\d3 



I1999I : iTorkelsson et aLlbOOd ). iLodato fc Pringlll|2007l )"find 
with smoothed particle hydrodynamics simulations of small 
and large-amplitude warps that the vertical viscosity satu- 
rates for small a such that 77^ < 3.5/a. Using these, the 
warp damping time scale becomes 



2r' ^ / Q \ ^r' 
U2 V3.5/ Ul 



2r' 



EC in5 a 115360 (^f'^'' 

5.6 X 10 — yr, (6) 



^3.5/ ui esF?7odd 

where 77odd = L/Lcdd- Equation ([5} can now be rearranged 
to give the radiation warping critical radius 

2rjl 2GA'h 



Ria.d > 



7c, 



(7) 



where 7crit — 0.32 (|Pringlel Il997l ). One way to obtain a 
stronger to rque is if the irradiatio n of the disc is driving 
an outflow jSchandl fc Meveill 19941 '). In this case, the torque 
is determined by the outfiow momentum or pressure at the 
sonic point, which is basically oc L, but also depends on 
the detailed disc structure in a complicated way. The en- 
hancement of the torque could be significant, as the ratio 
of momentum to energy for particles may be much larger 
than for photons. Although a gross over simplification, we 
may parameterize this by an additional (not constant) mul- 
tiplicative factor i^wind on the radiation torque, so that the 
warping criterion becomes 

^rad > 



277. 



^2 f 2 F2 
'crit wind 



2GAfbh „, in-3 -2 -2 „-2 

5 — = 9.1 X 10 a eo.i-Pwind PC, 



(8) 

where for the numerical value we have used 7crit — 0.32, 
e = O.leo.i and rj,^ ~ 3.5/a. Thus warping the disc in the 
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radial range of the GC discs (0.03-0.5pc) requires a ~ 0.3, 
or a ~ 0.1 and radiatively efficient accretion by a rotating 
black hole or modest enhancement of the torque by a disc 
wind. 



2.2.1 The Effect of Warping on Irradiation and 
Fragmentation 

Having argued for the possibility of radiation warping of the 
past accretion disc at the GC, we now discuss the feedback 
effect of warping on the irradiation and fragmentation of the 
disc.^ 

iPringld (|l997l ') studied the growth and non-linear evo- 
lution of radiation induced warps in the context of AGN 
accretion discs. His simulations show that highly warped 
discs exhibit limit cycle oscillations as a consequence of self- 
shadowing. As the warp gradually grows, the solid angle sub- 
tended by the central source decreases. When the maximum 
inclination of tt is reached, the instability cuts off leading to 
a decrease of the disc inclination. Thereafter, the disc ge- 
ometry again lets the instability set in, and the warp grows. 
Hence, even though the disc might be warped by large an- 
gles, radiation instability keeps acting on the disc so long 
as there is sufficient amount of radiation from the central 
parts. 

The e f fect o f warping on fragmentation is studied in 
iGoodmanI (|2003l ) who estimated for an irradiated disc the 
stabilizing temperature. Tat, above which fragmentation can 
not occur. When Tat is much larger than the temperature 
of the disc in equilibrium with radiation from the central 
source Toq, irradiation can not prevent fragmentation even 
for highly warped discs. The ratio Tst/Teq depends (weakly) 
on the mass of the central object and for a black hole of 
mass Mbh = 4 x 10^ Mq, and a disc with a warp of 60° it 
becomes ~ 1 at 0.03 pc and ~ 4 at 0.5 pc for a = 0.3, 
and e — 0.1. We note that for a disc around a black hole 
of mass Mbh = 10* M©, this ratio would have a value > 10 
at all radii, therefore fragmentation for a surrounding disc 
would be most likely. For the case of the GC on the other 
hand, the ratio Tst/Teq being close to unity at the inner edge 
of the disc suggests that fragmentation might be somewhat 
delayed until the effect of irradiation on the disc diminishes 
(see ^231l . 



2.3 Bardeen Petterson Effect and the Galactic 
Centre Stellar Disc 

In this section, we show under which conditions the Galac- 
tic Centre disc might have been warped due to Bardeen- 
Petterson effect. 

An accretion disc forming around a rotating (Kerr) 
black hole might initially have a total angular momentum 
misaligned with that of the black hole. Inner portions of 
the disc close to the black hole experience general relativis- 
tic frame dragging which causes a diffe rential precession, 
the s o-called Lense-Thirring precession (|Lense fc Thirrina 
Il918l ). As we have seen in the previous section, the viscous 
time scale increases with radius, hence the inner parts of 
the disc are forced to align with the black hole on time 
scales much shorter than those for the outer parts. Con- 
sequently, the disc develops a shape which for r < Rbp is 



aligned with the black hole, i.e. is flat, and for r Rbp 
its angular mo mentum direction changes g radually from ra- 
dius to radius (|Bardeen fc Petterson|[l975l ). Like the radia- 
tion pressure warping, the Bardeen- Petterson effect also is a 
competition betwee n the alignment time scale, a nd the pre- 
cession time scale (|Armitage fc Nataraianlll999l ). The pre- 
cession induced by the B ardeen-Petterson effect is given by 
iKumar fc Pringldll985l 'l 



(9) 



where a is the black hole spin parameter which can take 
values between and 1 for stationary, and maximally ro- 
tating black holes respectively. The precession time scale 
tbp = 2n/(j} is 



tbp 



(10) 



The critical radius, where the alignment time scale equals 
the precession time scale is then obtained by writing 

,,r2 = ^ -'^BP = ~ . (11) 

Using the steady state relations introduced in the previous 
section equation (|11[) becomes 



Rbp 



6aG^M^i,eE{RBp) 



(12) 



In order to estimate the Bardeen-Petterson radius for the 
GC, we need to solve equation p2|l for -Rbp- Writing 
T,{Rbp) = Eo/f"^ '* we obtain 



i?BP = 4.22 X 10 



-4 / aeQ/i536Q 

V ?;oddesF 



1/2.4 



pc. 



(13) 



This estimate shows that the Bardeen-Petterson radius for 
the GC is quite small, i.e. much below the inner edge of the 
observed discs. We should remind that the alignment time 
scale ti,2 for the assumed surface density profile is short. This 
means that the disc can be warped out to large distances 
since both Rbp and ti,^ depend on the disc parameters only 
weakly. 

The angular momentum of the disc at r is oc 
27rSrdrv^GMbh^ oc r^'^dr and varies only slowly with r. 
The ratio of the disc to black hole angular momenta can 
thus be estimated as 



Jd_ 

Jbh 



aGMlJc 



(14) 



where r^ — 2GMbh/c^ is the Schwarzschild radius, and 
Td ~ 0.1 pc. This suggests that by the time the disc is 
significantly warped at rd ~ 0.1 pc, the black hole spin 
should be significantly alig n ed: th e alignment time-scale 
found by iLodato fc Pringlei (|2006l ) under similar circum- 
stances is ~ 0.3ti/i ~ 0.3(3.5/a)J^2 ~ a~^ii/2- These argu- 
ments suggest that if the GC disc was first warped through 
the Bardeen-Petterson effect, and then fragmented to form 
the observed surface density of young stars, the star for- 
mation efficiency had to be high, esF < 1, to prevent the 
disc from dominating the angular momentum, and a < 0.1 
to give the disc time to warp before it is accreted and will 
align the black hole spin. 
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Figure 2. Normalized gravitational torque for the GC disc when 
the warp spans a range of —15° — > 15° in inelination. 



and gravity torques is 
Tr L 



^32g!?2^f!l, (19) 



^5360 



J 



6cr 27rGA/dE(r) J 

and the ratio of the viscous torque to the gravitational 
torque is 



-L =3.14 "^"""''^^ — . (201 

t grav Sec 27rGMdE(r)J ' Msaeo'^'^o.i J 

We see from Fig. [2] that the J-terms for the Galactic Cen- 
tre disc surface density profile are typically 0(0. 1). Thus at 
r = 0.1 pc Tr > Tg^av if esF < 1, but not if esF < 1- 
For a high central luminosity r^cdd and/or high star forma- 
tion efficiency esp the viscous torque on the disc dominates 
over its self-gravity. We have also seen for the case of the 
Bardeen-Petterson effect that esF < 1- This shows that it 
is justified to neglect the effects of the gravitational torques 
on the evolution of the disc in its active phase. 



2.4 Comparison of Gravitational, Viscous and 
Radiation Torques 



2.5 The Warped Stellar Disc After Fragmentation 



In ii2.2l we have discussed the conditions for radiation warp- 
ing. While equation ^ has to be satisfied for warping, it 
is also useful to compare the magnitude of radiation and 
gravitational torques acting on the disc. 

The magnitude of the radiation torque on a ring of 
radius r with radial w idth dr can be approximated as 
jOgilvie fc Dubusll200ll ) 



2TirdrTr — rdr, 

6cr 



(15) 



thus Tr — T,r^fl/tr- On the other hand the viscous torque 
on the ring, which tries to damp the warp in the case of 
radiation pressure warping, and which sets up the warped 
density distribution in the Bardeen-Petterson mechanism, is 
given by 

2nrdrn^ = ^1^1}^ = y.mrdr ^^^^ 



The gravitational torque of the disc on the same ring can be 
determined by integrating over the disc 



2nrdrT„ 



GMdT.(r)2irrdr 



J, 



(17) 



where J = T^mv / (GniiMd/ri) is a dimensionless integral 
depending on r/ri, r/r2 

27vr'dr'Eir') r\' , 9/3 

^ = X^ Md ir-+r-f/^ s.nmi{P,r/r)-, 

(18) 

and where /3(r, r') is the angle between the normals of the 
two rings at r and r', is the inclination of the ring at r, 
so dp /d o = 0(1), and 1(13, r/r') is the integral in equation 
(13a) of lUlubav-Siddiki et al.l (|2009l ) (hereafter US09). 

We show in Fig. [2] the values of the J-terms at different 
distances from the black hole when the warp spans a range 
of —15° — ^ 15° in inclination. For larger warps, the gravi- 
tational torques are weaker. Also if we include the part of 
the disc outside [rin, rout] the product M^J decreases as the 
torques are dominated by the nearby parts of the disc. Using 
the surface density profile above, the ratio of the radiation 



In sections (|2.2p and (|2.3p we have seen that for a range of 
assumed accretion disc parameters, the disc at the GC could 
have been warped in the accretion phase. The evolution of 
the disc is likely to be governed by the viscous or radiation 
torques during this phase. 

We now assume that some time after the disc has be- 
come warped, the active phase of the GC ends. Since ?7cdd 
and the accretion rate will then be highly reduced the disc 
will receive much smaller energy input and it is reasonable 
to assume that it can now cool rapidly. If the cooling time 
is short, fragmentatio n and star form ation can occur on a 
dynamical time scale (|Gammie|[200lh . The gas content of 
the disc could then be converted into stars with efficiency 
esF and the rest of the gas will be lost. 

Subsequently, the now stellar disc would be subject only 
to gravitational forces, and precess accordingly. In addition 
to the precession caused by the disc's self-gravity, the sur- 
rounding old star cluster will also contribute unless it is 
spherically symmetric. After analyz i ng th e proper motions 
of the old cluster stars, iTrippe et al.l (|2008l ) reported a small 
amount of rotation in the direction parallel to the Galactic 
ro tation. The ro t ation of the old cluster is also confirmed 
bv lSchodel et^ (|2009l l. Based on these findings, one might 
consider the possibility of the star cluster being slightly fiat- 
tened, hence have an effect on the long term evolution of the 
disc. 



3 WARP MODEL FOR THE GALACTIC 
CENTRE DISC 

In the previous sections we have argued that for a range of 
parameters, a flat disc around SgrA* could become warped 
before forming stars, after which the disc evolves gravita- 
tionally. In this section, we will introduce the model we have 
used to follow the time evolution of the warped stellar disc. 



3.1 The Equations of Motion 

As in US09, we model a warped disc as a collection of con- 
centric circular rings which are tilted with respect to each 
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other. The rings are in gravitational interaction with each 
other, and are also torqued by the surrounding old star clus- 
ter. They are characterized by their masses rrn, and radii n. 
The fast orbital motion around the SMBH of mass Mbh is 
thereby time-averaged. The geometry of the rings is defined 
by the Euler angles {ip,8,(j>). For a system of n rings, the 
equations of motion for any of the rings i are given by 



rriirf 



(Pi sm di cos t 



iP^ . sin Qi 



dVi 



d(t)i 



(21) 



(22) 



(23) 



(24) 



Here p^. = mirlQ,i is the orbital angular momentum, (^i 
and 9i are the rates of precession and nutation caused by 
the torques dVi/dOi and dVi/d(j>i respectively, where Vi is 
the potential energy of ring i in the field of the other rings, 
and of the old stellar cluster. The details of the calculation 
of the mutual torques can be found in US09. To calculate 
the torque due to t he surrounding star cluster, we use the 
description given in ISparkd (|l986l ) appropriate for slightly 
flattened systems. The potential energy of a ring in the field 
of the cluster is given by 



Vci = —^(1 - q) sm ( 



(25) 



where vq is the assumed constant circular velocity generated 
by the star cluster, and g is the fiattening. The torque on 
ring i due to the star cluster is then written as 



rriiVo 



^ = ^(l-g)sm2e. 



(26) 



The precession frequency of the disc induced by the star 
cluster at different radii r is given by 



VQ 

Sri 



(1 — q) cos ( 



(27) 



3.2 Numerical Setup 



The numerical approach we use to study the time evolution 
of the disc is similar to the one presented in US09. There, we 
integrated the equations of motion for the rings after evalu- 
ating the steadily precessing equilibria for the given param- 
eters to test the stability of those configurations. For the pa- 
rameters relevant to the GC disc we do not find such steadily 
precessing equilibria. In an earlier work it is shown that cen- 
trally illuminated discs around black holes might under cer- 
tain ci rcumstances b e warped by large angles, e.g. even by 
~ TT (iPringldl 19971 ) . The final geometry of the warped disc 
depends on the model assumptions like the strength of the 
perturbation leading to radiation instability and the loca- 
tion in the disc where the perturbation occurs. For the case 
of Bardeen-Petterson warping on the other hand, the de- 
gree of the warp depends on the initial mutual inclination 
between the angular momenta of the gaseous disc and the 
black hole. Therefore the most realistic initial disc shape for 




0.2 0.3 
radius (pc) 

Figure 3. Initial inclination profile of the model discs. 



our simulations can only be chosen after carrying out sim- 
ulations of gaseous discs subject to radiation instability, or 
Bardeen-Petterson effect, which is beyond the scope of this 
study. Therefore, without attempting to find such precise 
initial conditions, we adopt as initial warp shapes, the disc 
configurations which are likely to be imposed by the warping 
mechanisms discussed in the previous sections. Specifically, 
all the models are assigned an initial warp profile such that 
the inclinations range from —30° to 30°, and the azimuthal 
angles are zero. We note that this choice of the inclinations 
is in the range of the so-far reported values for either of the 
mechanisms. Fig. [3] depicts the initial inclination profile for 
the models we constructed. 

In our simulations, we model a disc as a collection of 
50 equally spaced, circular rings. The model discs extend 
between 0.03 — 0.5 pc, surrounding a black hole of mass 
A'/bh = 4 X l(fMQ. The surface density of the discs de- 
cline as l/r^ * as indicated by the observations. Initially, 
the discs are given a precession frequency such that when 
normalized to the orbital frequency of the innermost ring 
it gives 4>/flin ~ —10^^. Below we will often give evolu- 
tion times in terms of the orbital time at the inner edge 
of the disc, (Porb)in ~ 244 yr. We consider discs which are 
under the influence of only self-gravity (models SG), and 
discs which experience the effects of the star cluster and the 
self-gravity together (models GLSG). We assume three dif- 
ferent mass fractions, Md/Mbh =0.00134, 0.0134, and 0.134 
(LM-low mass, IM-intermediate mass, and HM-high mass 
respectively). We note that assuming Md/Mbh = 0.00134 
implies a star formation efficiency of esF = 1, a value for 
which warping the planar GC disc is most plausible (see 
321 . However, in order to keep our discussion broad, and to 
make our models applicable to systems other than the GC, 
we also explore mass fractions which are higher than that 
inferred from the observations of the GC discs. 

Observations of the nuclear star cluster of the GC sug- 
gest a mass range of M* ~ [1 — 2] x IO^Mq w i thin the outer 
edge of the discs at 0.5 pc (|Trippe et al.|[200^ : ISchodel et al.l 
l2009i ). Therefore, for the GLSG models, we set the quantity 
vq{1 — q) that determines the strength of the cluster torque 
to 860 km^ s~^ (see equation (|26[) V This value can be ob- 
tained, for example by setting the cluster mass within 0.5 
pc to lO^Af©, hence adopting a circular velocity of no ~ 93 
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(106 yr) 


Break Up 


SG.HM 


0.13400 


0.03-0.5 


60 
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0.6 


yes 


SGJM 


0.01340 


0.03-0.5 


60 


10-5 






8 


yes 


SG.LM 


0.00134 


0.03-0.5 


60 


10-5 






60 


no 


CLSG.HM 


0.13400 


0.03-0.5 


60 


10-5 


860 




0.35 


yes 


CLSGJM 


0.01340 


0.03-0.5 


60 


10-5 


860 




0.8 


yes 


CLSG.LM 


0.00134 


0.03-0.5 


60 


10-5 


860 




1.6 


no 





















Table 1. Parameters of the simulations discussed in the paper. The first column shows the label of the simulation. SG models are those 
considering only the self-gravity of the disc, while CLSG models include the effects of the star cluster as well. HM, IM, and LM stand for 
high mass, intermediate mass, and low mass respectively, which refer to the mass fractions given in column two. The third, the fourth, 
and the fifth columns show the adopted radial range, the overall warp, and the normalized initial precession frequency respectively. The 
sixth column gives the value of the quantity Vq{1 — q) which measures the strength of the torque from the star cluster. The seventh 
column lists the precession times attained at the end of each simulation, where the last column summarizes the behavior of the disc in 
terms of its coherence. The underlined parameter values follow from the observations or the models, whereas the others are assumed 
values. 



km/s as induced by the star cluster, and the flattening of 
the cluster to q — 0.9. 

The parameters of the simulations along with the result- 
ing precession time scales and evolution are listed in Table [T] 
In these simulations, the total energy of the rings is typically 
conserved to an accuracy of 10-5. 



4 RESULTS 

In this Section, we describe the main results that emerged 
from our simulations. We discuss in turn the precession of 
the model discs and the evolution of the discs' warp shapes. 



4.1 Precession of the Disc 

The self-gravity and the star cluster torques will force the 
disc to precess on time-scales proportional to the disc mass 
fraction, and the quantity Wo(l— g) respectively (see equation 
(28) in US09 for the case of purely self-gravitating discs). 
The overall evolution of the disc will then be determined by 
the relative strengths of both torques. 

In Fig. [4] we show for the SG models how the preces- 
sion frequencies evolve during the simulations. The x-axis 
shows the elapsed time in terms of the orbital frequency of 
the innermost ring on the lower x-axis, and in units of 10^ 
yr on the upper x-axis. The y-axis shows the absolute value 

n 

of the averaged precession frequency, i.e. |^^(^i/n|. The 

red (upper), the black (middle) and the blue (lower) curves 
correspond to disc mass fractions of Afd/A/bh = 0.13400, 
Md/Mbh = 0.01340, and Afd/Afbh = 0.00134 respectively. 
We see that although the average precession frequencies ex- 
hibit small oscillations in amplitude, they attain their true 
values (imposed by the disc's gravity) soon after the simula- 
tion starts, remaining approximately constant even though 
the shape of the disc evolves, as is shown in the next sub- 
section. While individual rings may precess in the forward 
direction for short time intervals, most rings precess back- 
wards for most of the time, so that the average precession is 
always retrograde. 

Throughout the simulation model SGJHM precesses 
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0.5 1 1.5 2 



2.5 



10' 
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Figure 4. Evolution of the precession frequencies of the models 
considering only self-gravity of the disc. Different curves corre- 
spond to models with different mass fractions. The precession 
frequency increases in proportion to the disc mass fraction. 



with a normalized frequency of about — 4x 10^'', correspond- 
ing to a precession time scale of Tp ~ 6 x 105 yr. For model 
SGJM the normalized precession frequency drops to approx- 
imately —3 X 10^5 giving Tp ~ 8 x 10^ yr, and when the mass 
fraction is decreased even further we see from model SG_LM 
that 0/r2in = — 4 X lO-*" leading to Tp « 6 x lO'^ yr. (In or- 
der to check this, we have computed for model SG_LM the 
normalized precession frequencies that the rings would have 
under the self-gravity torques when the initial warp profile 
is considered (see equation (14) in US09). The values are 
in the range [—8,-2] x 10-^, with some of the rings hav- 



ing (j)/ilin 



-4 X 10 , consistent with Fig. g]) We can 



readily conclude from this, that the precession frequencies 
scale with disc mass. In the equations of motion, when the 
quadratic term for the precession frequency is set to zero, 
all terms scale linearly with ring mass so the precession fre- 
quency should indeed be proportional to the mass. Among 
the models discussed above only SGJHM would have com- 
pleted a full precession during the life time of the stars. 

In Fig. [5] we display how the above scaling manifests 
itself when we look at the averaged azimuthal angles of the 
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time (Porb)ii 



Figure 5. Precession of tlie models considering only self-gravity 
of the disc. Different curves correspond to the models with dif- 
ferent mass fractions: the red (upper) curve is for model SG_HM, 
the black (middle) curve is for model SG_IM, and the blue (lower) 
curve is for model SG_LM. The ever faster precession of the more 
massive discs can be clearly seen. 



rings, where the averaging is done in a similar way as for 
the frequencies. The different curves correspond to the dif- 
ferent mass fractions again. At the end of the simulation, 
10000 orbital periods after the simulation has started, model 
SG_HM is seen to have precessed approximately 4.5 times, 
while model SGJM has completed nearly half a precession. 
The precession of the least massive model SGJLM on the 
other hand can barely be noticed with 0.05 times a full pre- 
cession. 

We now move on to the models including the effect of 
a surrounding cluster. In Fig. [5] we show the evolution of 
the averaged precession frequencies for the CLSG models. 
In contrast to the SG models, the precession frequencies for 
the CLSG models may not reach a constant value (note that 
the noise in the curves is the same as in Fig. [4l given the 
axis scaling). Until 3000 orbital periods, the frequencies of 
all the models increase linearly in time. Since the param- 
eters of the star cluster are the same for all the models, 
it is again the most massive disc which precesses with the 
highest rate 4>/Q.in = —7 x 10~* corresponding to a pre- 
cession time scale Tp ~ 3.5 x 10^ yr. Comparing models 
SG_HM and CLSG_HM, we see that the averaged preces- 
sion rates are approximately the same for the first 1000 or- 
bital periods. At the end of the simulations, the star cluster 
leads to an approximately 2 times faster precession than in 
the self-gravitating case when M^/My^i-^ = 0.134. For model 
CLSGJM the precession is slower due to the decreased disc 
mass, but attains a value of 4>/Qin ~ —3 x 10~* giving 
Tp Ri 8 X 10^ yr, 10 times shorter than for model SGJM 
due to the influence of the star cluster. Finally the lowest 
mass disc model including the star cluster CLSG_LM has 
<})/0,in = —1.5 x 10~* and Tp ~ 1.6 x 10^ yr, which is nearly 
40 times shorter than for SG_LM. 

It emerges from these comparisons that the effect of the 
star cluster on the disc is most pronounced when the disc 
mass is low, otherwise the self-gravity torques dominate the 
evolution of the disc. In Fig. [7] we show for different mass 
fractions the expected ratio (based on equation (|26|) and 
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Figure 6. Evolution of the precession frequencies of the models 
including the torques from self-gravity and from the surrounding 
star cluster. Different curves correspond to models with different 
mass fractions as indicated on the labels. 
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0.1 0.2 0.3 0.4 0.5 

radius (pc) 

Figure 7. The comparison of the cluster torques, Td, to the self- 
gravity torques Tgrav acting on the rings. The low mass disc is 
dominated by the torques from the star cluster, while the high 
mass disc is dominated by self-gravity torques. 

equations (13) in US09) of the cluster torque, T^, to the self- 
gravity torque, Td, at different locations in the disc when 
the initial warp profile is considered. One can see that the 
low mass disc is indeed dominated by the cluster, while the 
high mass disc is dominated by self-gravity except in the 
inner parts. The kink and drop around ~ 0.25 pc is caused 
by the change of sign of the torques which follows from the 
choice of the warp profile. 

The relative effects of the cluster and the self-gravity 
torques can also be observed in the azimuthal profiles, i.e. 
radius versus azimuth, of the CLSG models which we show 
in Fig. [8] The black horizontal dotted line shows the ini- 
tial profile of the azimuthal angles. The red curve with filled 
circles shows the azimuthal profile of model CLSG_IIM to- 
wards the end of the simulation. The black curve with filled 
triangles, and the blue curve with filled squares depict the 
profiles of models CLSGJM, and CLSG XM respectively. 

In general, it is possible to find disc configurations where 
the azimuthal angles of the rings are aligned even in the pres- 
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Figure 8. Azimuthal profiles of the CLSG models showing pre- 
cession angles for all rings after 6000 inner periods. The black hor- 
izontal dotted line at the top shows the initial configuration. The 
red curve with filled circles shows the azimuthal profile of model 
CLSG_HM. The black curve with filled triangles, and the blue 
curve with filled squares depict the profiles of models CLSG_IM, 
and CLSG_LM respectively. The effect of the star cluster on the 
disc is best observed for model CLSG_LM where the azimuthal 
angles exhibit a differential profile. 

ence of a surrounding star cluster (US09), however this can 
only be achieved for certain combinations of the disc param- 
eters. Although for the GC parameters such configurations 
were not found, we see from Fig.[5]that an almost steady pre- 
cession can be achieved when the disc is massive enough as in 
model CLSGJHM. In this case, the inner and the outer parts 
of the disc are locked together, while the middle parts which 
are ahead of the others are also seen to precess together. For 
model CLSGJM where the effect of the self-gravity torques 
is diminished, the star cluster torques start to dominate in 
the inner parts of the disc, where the azimuthal profiles show 
more of a differential profile rather than a flat one. Since the 
rate of precession induced by the star cluster decreases go- 
ing outwards in the disc (see equation (|27|l '). the outer parts 
for this mass fraction are still held by the self-gravity of the 
disc. For the lowest mass model CLSG_LM, we see that the 
star cluster torques dominate over the self-gravity torques 
and the rings making up the disc can not precess with a 
common azimuth. 

4.2 Evolution of the Inclinations and Breaking Up 
of the Disc 

In this section we will discuss how the initial warp shape 
of the disc is modified under the self-gravity and the star 
cluster torques. Again, if we had started the simulations 
adopting equilibrium inclinations, we could expect to see 
that for a range of parameters, the discs would be stable, 
hence would not change their inclinations. Since we have 
seen in the previous section that even in the absence of an 
equilibrium configuration, some discs can maintain a near 
steady precession -at least at the inner and outer parts of 
the disc-, we might expect to see some behavior mimicking 
equilibrium in the the long term evolution of the inclinations 
as well. 

In the plots we will be depicting in this section, the 
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Figure 9. Time evolution of the inclinations for model SG_HM. 
The disc breaks into separate pieces which become pronounced 
after 330 orbital periods (marked by the vertical red (left) line). 
After 890 orbital periods (marked by the vertical blue (right) line) 
the inclinations reach a maximum separation of about 40° . 



X-axis will show the elapsed time in terms of the orbital 
period of the innermost ring on the lower x-axis, and in 
units of 10^ yr on the upper x-axis. The y-axis will show the 
ring inclinations, where the negative ones correspond to the 
inner parts, and the positive ones correspond to the outer 
parts. At any instant of time t, a vertical cut through the 
plot will then give the warp profile, i.e. the inclination at all 
radii. 

We start our discussion by considering first the SG mod- 
els. In Fig. [9] we show the evolution of the individual ring 
inclinations for model SGJIM. We see that soon after the 
simulation starts the disc breaks into pieces, where the rings 
lying in the middle parts are pushed to more polar orbits, 
and the edge rings become more equatorial. After 330 or- 
bital periods (indicated by the red (left) vertical line), the 
breaking of the disc becomes very prominent, leading to a 
configuration of three concentric, separate discs, the middle 
one being the less warped one. At about 890 orbital pe- 
riods (indicated by the blue (right) vertical line), the two 
positively inclined discs merge together, and the separation 
between the inner and the outer discs reaches a maximum of 
about 40° . At this phase, the inner disc becomes almost pla- 
nar, while the outer one is warped by about 10°. At about 
1480 orbital periods the inclinations of the inner rings ap- 
proach their initial values, and the mutual inclination of the 
separated discs become ~ 25°. After this time, subsequent 
oscillations with almost equal amplitude and period follow. 
Although a more detailed investigation of this behavior is 
beyond the scope of this study, we should note that these os- 
cillations are reminiscent of stable oscillations for which the 
system parameters oscillate around their equilibrium points 
once perturbed. 

In Fig. [To] we show the warp evolution of model SGJM. 
We see again that the disc breaks into pieces. For this model 
however the breaking becomes pronounced at about 3300 
orbital periods (indicated by the red (left) vertical line), 10 
times longer than for model CLSG_HM corresponding to the 
10 times longer precession time. At this stage, the disc has 
broken into three pieces, with similar characteristics, i.e. the 
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Figure 10. Time evolution of tiie inclinations for model SG_IM. 
The evolution is very similar, but retarded by a factor of 10, writli 
respect to the evolution of model SG_HM. 
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Figure 11. Time evolution of the inclinations for model SG_LM. 
The inclinations stay almost unaltered during the simulation 
which is an effect of the low mass fraction adopted. 



degree of the warps, and the separation between the other 
discs, to those in model CLSG_HM. At about 8900 orbital 
periods (indicated by the blue (right) vertical line) , the two 
of the positively inclined discs combine to form a single one, 
which is again separated by 40° from the negatively inclined 
one. 

The warp evolution of model SGJLM is shown in Fig. 
111! The inclinations of the rings for this model are seen 
to remain almost unchanged during the simulation. This is 
an effect of the low mass adopted for this model, and the 
corresponding long precession time. In order for this model 
to show a prominent breaking, the disc would have to evolve 
for about 3.3 x 10* orbital periods (~ 8 x 10^ yr), which is 
longer than the age of the stars. 

In Fig. [12] we show the warp evolution of model 
CLSG_HM which considers the effect of both the self-gravity 
and the star cluster. We have seen in the previous subsection 
that for this mass fraction, it is the self-gravity that mainly 
dominates the precession evolution of the disc. The situation 
is therefore similar to that in Fig. [9] when we consider the 
warp evolution. We see that model CLSG_HM also breaks 
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Figure 12. Time evolution of the inclinations for model 
CLSG_HM. The disc breaks into pieces on time scales comparable 
to those for model SG_HM. The long time behavior of these two 
models differ in that when the cluster is present, the oscillations 
in inclination are less well behaved. 



into pieces on time scales comparable to those for model 
SG_HM, attaining similar separation angles. More specif- 
ically, model CLSG_HM reaches a maximum separation of 
40° at about 600 orbital periods. At 1300 orbital periods, the 
mutual inclination between the broken parts drops to 10°. 
This is twice the value attained in model SG_HM. However, 
the major difference between the two models shows itself in 
the long term evolution, the cluster model exhibiting a less 
well behaved oscillation pattern. 

When the mass fraction is decreased, the effect of the 
cluster torques become more prominent. This can be seen in 
Fig. 1131 After about 500 orbital periods the ring inclinations 
start to deviate from their original values. 

Compared to the other models in which the disc breaks 
up, the evolution for this model is more disordered. There is 
no well defined pattern left in the warp profile once the disc 
breaks. After 2500 orbital periods, the positively inclined 
rings are pushed to inclinations of more than 30°, reaching 
a maximum of about 50° at the end of the simulation. At 
this stage, only some of the inner rings exhibit a disc-like 
structure with a warp of about 10°. The original disc as a 
whole has dissolved due to the differential precession induced 
by the cluster, and the self-gravity of the disc is not able to 
hold the disc intact. 

In Fig. [14] we show the evolution of inclinations for 
model CLSGJLM. For this model, we see that the inclina- 
tions are altered by only a few degrees, which at the end 
of the simulation leave the overall warp shape almost un- 
changed. This is because for this model the torques are dom- 
inated by the star cluster, leading to nearly free differential 
precession of the rings at constant inclination, causing the 
rings to spread in azimuth by about one full rotation (see 
Fig.[H]). 

In the last part of this section, we show in Fig.llSlseveral 
examples of the 3-dimensional visualizations of our models 
at different stages of their evolution. The color code displays 
the height (z-coordinate) of the disc. Different viewing an- 
gles are chosen to highlight the various structures present in 
the warp shapes. The figures correspond to: the initial con- 
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After application of proper selection criteria, the sam- 
ple includes a total of 136 stars on clockwise and counter- 
clockwise rotating orbits. These stars have measured {x,y) 
position, while the line-of-sight distance z is unknown. For 
the analysis, z is therefore treated with a Monte-Carlo sim- 
ulation. The respective space ve locities Vx, Vy, and Vz are 
also given in iBartko et all (|2009l '). 

The n umber density of st ars in the discs, Enr(''), scales 
as (|Bartko et al.llioiol '). Here r is the 3-dimensional 

distance to SgrA*, corresponding to the radius r of our rings. 
The total number of stars on a ring is then 27rrArEnr(r). If 
we assume that the total number of stars in the disc is A*'t 
we can write 



^27rrArEnr(r) = iVt, 



(28) 



Figure 13. Time evolution of the inclinations for model 
CLSGJM. The self-gravity of the disc can not hold the rings 
together, and the disc breaks up into several fragments under the 
influence of the star cluster. 
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Figure 14. Time evolution of the inclinations for model 
CLSG_LM. The inclinations do not change substantially while 
the rings precess apart azimuthally (see Fig.[8ll, because the mass 
fraction adopted is low for this model. 



figuration common to all the models (a), model SGJHM at 
t = 340 (b), model SG.IM at t = 8900 (c), model CLSGJM 
at t = 2500 (d), model CLSGJM at t = 5680 (e), and model 
CLSG_LM at t = 5680 (f) respectively where t = t/Porbj,,. 
The broken shapes of the discs after some evolution for mod- 
els SGJIM and SGJM are apparent in figures (15b) and 
(15c). Fig. (15e) depicts clearly how the shape of a disrupted 
disc with parameters identical to those for model CLSGJM 
would look like. In the last panel in Fig. (15f), we see that 
the inner parts of the disc precessed more compared to the 
outer parts, while the inclination of the disc has only slightly 
changed compared to its initial state (see figures (O and 



d) 



4.3 Comparison With the Observations 

In this section, we compare our time evolution models with 
the observations of the GC discs. In doing so. we will make 
use of the results of iBartko et~all (|2009l . I2OI0I '). 



where n is the number of rings making up the disc. We 
write Enr('') = Eo(r/rin)~^ *. The normalization Eo, i.e the 
number of stars on the innermost ring, is then obtained by 
writing 



En = 



Nt 



1 

27rAr ri^* V 



(29) 



The coordinates of a star on a ring are given by 

X = r{cos cj) cos ill — cos 6 sin (f) sin ip), 
y = r(sin (/)cos ?/) -I- cos 5 cos (;/>sin 7/)), 
z = r sin S sin 7/). 



(30) 



In order to construct the phase space distributions of 
the stars populating the discs, we write the components of 
the stellar velocities 

Vx ~ r{0 sin 9 sin il^ sin (f) — tp sin ip COS (f) — (p sin (j) COS 
— cos ^ cos sin t/i — t/j cos S sin (;/>cos 7/)), 

Vy — r{(f>cos<l)costp — (j)cosOsin(j)sintp 

+ xj) cos 6 cos (j) cos tp — 6 sin 6 sin ^ cos (j> — ip sin tp sin <j)) , 

Vz — r(^ sini/; cos S 1/) sin S cos i/)). (31) 

Since we know the values of (0, 9, tp, (p, 6) from our time inte- 
gration, we can generate a random sample of stars on these 
rings, by dicing ^. We can then construct the phase space 
distribution of stars at each instant of time t using equations 
(|30ll and pTI . 

The GC contains two warped and almost orthogonal 
discs of youn g stars which extend over a similar range of pro- 
jected radii (jBartko et al.ll20l"ol l. The CW disc is observed 
to extend between projected radii of approximately 0."8 and 
15", while the CCW disc is seen between projected radii of 
about 3." 5 and 15". This makes it unlikely that the two 
stellar discs derive from a single progenitor disc, but rather 
suggests that they might be tracers of two separate accre- 
tion episodes. In the following, we will concentrate in our 
model comparison on the CW disc, which is better defined 
due to its larger number of constituent stars. 

We have shown in the previous section that massive 
discs break into pieces. Although the broken discs may 
precess as coherent separate discs under certain circum- 
stances, their degree of warping is highly reduced compared 
to the original disc. The observational analysis presented in 
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Figure 15. 3-dimensional visualizations of the warped disc models at different stages of their evolution. The color code displays the 
height (z-coordinate) of the disc. Different viewing angles are chosen to highlight the various structures present in the warp shapes. The 
figures correspond to: the initial configuration common to all the models (a), model SG_HM at t = 340 (b), model SGJM at t = 8900 (c), 
model CLSGJM at t = 2500 (d), model CLSG.1M at t = 5680 (e), and model CLSG.LM at t = 5680 (f) respectively where t = t/Porbi„- 



iBartko et al] (|2009l ) can discriminate between a disc precess- 
ing as a single unit and a broken disc, such that the former 
produces a warp profile that changes gradually between the 
inner and the outer radii, while the latter would result in a 
jump at a certain radius. 

A useful way of comparing the models with the obser- 
vations is to construct the components of angular momenta 
and see how they change with projected radius. The compo- 
nent of angular momentum along the axis pointing towards 
the observer will then give information on the twist of the 
disc, while the component projected onto the plane of the 
sky will give the degree of the warp (for further de tails on 
the definitions of the angles see iBartko et all (|2009l )'). 

Fig. Uni shows the component of the local average angu- 
lar momentum direction of the m odel discs (red tria ngles), 
and of the CW disc as published in lBartko et al.l (|2009l ') (blue 
dots) on the plane of the sky as a function of average pro- 
jected distan c e from the centre {\—(j}j{R) in the notation of 
IBartko et al.l (|2009l ) where J denotes the angular momen- 
tum per mass). The data from the simulations is rotated 
successively such that the orientations of the model and the 
observed discs match at the inner edge, and the warp in the 
simulations fit the observed warp best. The points referring 
to the simulations are chosen at certain times during the 
late stages of the evolution of the discs. The figures on the 
top are for the self-gravity-only models, and the ones at the 
bottom are for the models including the effects of the star 
cluster. The mass fraction of the simulated discs decreases 
going from left to right, therefore the models should be read 
as HM, IM, and LM going in this direction. 

We see that a good agreement with the observations 



of the CW disc is achieved for our model SGJLM in that 
the inclination (warping of the disc) changes gradually be- 
tween the inner and the outer edges of the disc, and that 
the amplitude of the warp is successfully reproduced. For 
the self-gravity-only models, higher mass fractions lead to 
flattening of the inner parts of the disc, hence can be ex- 
cluded to be representative of the data. On the other hand, 
when the torques from the star cluster are taken into ac- 
count, none of these models seem to match the data well. 



5 DISCUSSION 

5.1 Dependence on the Initial Conditions 

In the previous section we argued that a good agreement 
with the observations of the CW disc is achieved for our 
model SGXM. To test the importance of the choice of the 
initial warp configuration we performed a run, model LI, 
where the initial inclinations, i.e. A6, span a range —15° to 
15° , and the other set of parameters are identical to those of 
model SG_LM. Fig. 1 171 shows the evolution of the azimuthal 
angles for models LI (solid line), and SG_LM (dashed line). 
We see that when the disc is less inclined it precesses faster, 
by a factor of about 1.3 for the parameters chosen. An an- 
alytic expression showing the inverse proportionality of the 
precession rate to the disc inclination for steadily precessing 
discs can be found in US09. 

We show in Fig.[TS]the evolution of the ring inclinations 
for model LI. We see that the splitting of the disc into parts, 
although not much pronounced, occurs on a shorter time 
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Figure 16. The warp profiles of the model discs (red triangles) and of the CW disc (blue dots). The figures on the top show the 
self-gravity-only models, and the ones at the bottom show the models including the effects of the star cluster. The mass fraction of the 
simulated discs decreases going from left to right. The lowest mass model for the case of self-gravity-only gives the best agreement with 
the data. 
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Figure 17. Time evolution of the azimuthal angles for models 
LI (solid line), and SG_LM (dashed lino). Model LI precesses 1.3 
times faster than model SG_LM due to its decreased overall in- 
clination. 

as compared to model SGXM (the inclinations in model 
SGJLM had remained unaltered during a similar time span). 
Even if there are slight excursions of the ring inclinations 
from their initial values, the overall warp shape does not 
change considerably throughout the simulation for model 
LI either. Hence we conclude that in order for the low mass 
disc to evolve towards the observed orbital configuration, its 
initial inclination determined by the warping mechanisms 
has to be close to the observed inclination. 



Figure 18. Time evolution of the ring inclinations for model LI. 
The ring inclinations change only slightly throughout the simula- 
tion, leading to a final warp of about 32° . 



5.2 Comparison to Linear Evolution 

In US09 it is shown that linear theory over estimates the 
torques acting on the edges of the disc, leading to a slightly 
less inclined warp configuration at its maximum possible 
value (within the validity of the linear theory). Here, we 
would like to briefly discuss the evolution of the disc in lin- 
ear theory, where the linearization is done as in US09. In 
order to reduce computational time we have considered one 
purely self-gravitating model adopting a mass fraction of 
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Figure 19. Average precession frequencies attained during the 
simulations for a mass fraction of Mj/Mbh = 0.134. The sohd line 
is for model SG_HM calculated in the non-linear torques regime, 
and the dashed line is for the linear model. The disc processes 
with a much higher rate in the linear regime. 

A^d/A^bh ~ 0.134. We have run the simulation for about 
1250Porbi„ 1 which corresponds to half a precession time for 
the non-linear model for this mass fraction. Since the evo- 
lution of the disc scales approximately linearly in mass, the 
comparison which we will present should equally apply to 
discs with other mass fractions. 

In Fig. [1^] we show the averaged precession frequencies 
of the non-linear (model SG_HM) (solid line), and the linear 
models (dashed line). We see that in linear theory the model 
acquires an average precession frequency of 4>/flin ^ — 6 x 
10~^, corresponding to a precession time scale of 4 x 10'' 
yr. Hence the model disc in linear theory precesses 15 times 
faster than it should. Turned around, the system behaves 
as if it were 15 times more massive (for the mass fraction 
considered here, this would unphysically indicate a disc twice 
more massive than the central black hole). 



5.3 Limitations 

5.3.1 Eccentricity of the Disc 

In our simulations, we used a simple circular ring model al- 
though the observations point to an eccentricity of the orbits 
with e ~ 0.35 for the C W disc and close to zero eccentric- 
ity for the CCW disc (jPaumard et al.l l2006l : iBartko et al] 
I2OO9I . I2OIOI ). The change of eccentricity (as well as incli- 
nation) in an initially circular disc might be driven by 
the interaction of the disc stars with the s urrounding cusp 
llPerets et al. I I2OO8I : iLockmann et al.1 120091 '). Following the 
time evolution of single and ini t ially mutually inclined two 
discs. [L5ckmann fc Baumgardtl (|2009l ) showed that such in- 
teractions can cause the inner parts of the disc, i.e. r < 0.3 
pc, to gain quite high eccentricities on a typical time scale 
of 5 Myr. Neg l ecting the effects of the disc's self-gravity 
iMadigan et al T (|2009t ) showed that an initially eccentric disc 
embedded in a stellar cusp around a massive black hole is 
subject to an instability which leads to the growth of ec- 
centricities on about a precession time (~ 0.6-0.7 Myr at 
0.05 pc). It is not known whether taking into account the 



self-gravity of the disc would modify the results of these 
calculations strongly. 



5.3.2 Viscous Evolution of The Disc 

In proposing our scenario, we made a strong assumption 
that there had been an accretion disc around SgrA* a few 
million years ago. While we have argued for the plausibility 
of disc warping prior to star formation, our simulations were 
performed following only the gravitational evolution of the 
disc. The next step forward would be to consider the vis- 
cous evolution of the disc under the non-linear self-gravity 
torques, as well as the torques from the old star cluster. As 
we will briefly discuss in the next subsection, viscous discs 
in various physical environments are also prone to be broken 
or disrupted by differential precession. 



5.4 Vector Resonant Relaxation as a Warping 
Mechanism 

In our scenario, we considered two mechanisms, the radi- 
ation instability, and the Bardeen-Petterson effect to ex- 
plain the warped gaseous disc origin of the GC discs. With- 
out making a detailed numerical analysis of these mech- 
anisms, we argued for their possibility based on simple 
arguments such as the critical radii and time scales for 
warping. It would be interesting to explore other warp- 
ing mechanisms than the t wo discussed here. Recently, 
iBregman fc Alexander] (|2009l ) showed that vector resonant 
relaxation (VRR) resulting from the interaction of the disc 
with the surrounding old star cluster might explain the warp 
in the maser disc of NGC4258. When scaled to the param- 
eters of the GC region, VRR serves as a viable mechanism 
for o rbits within the central ~ 0.2 pc (|Kocsis fc Tremaind 
Beyond this radius, the VRR time scale exceeds the 
life time of the disc stars. However, it would still be inter- 
esting to see how the disc responds to the stellar cusp when 
the simulations are carried out in the non-linear regime. 



5.5 Broken/Disrupted Warped Discs 

We have shown examples of purely gravitating warped disc 
simulations where the discs with a certain mass fraction 
break into pieces. This is an effect of the faster precession of 
some parts of the disc, and can lead to a complete disruption 
of the disc if the broken parts can not be held together by 
internal torques. 

Broken or disrupted discs have been reported before 
mainly in the context of viscous discs. Accretion discs 
warped by tidal effects in misaligned binaries can be dis- 
rupted due to differential precession when the disc as- 
pect ratio h = H/r ~ 0.01, where H is the disc height 
(|Larwood et al.lll996l : iFragner fc NelsonI [2OI0I ). When the 
warp propagation through the disc is diffusive rather than 
wave-like, isothermal discs with small viscosity exhibit a 
steepening in their warp profiles where the dis c breaks into 
two parts (|Ogilviell200"g : iLodato fc Pricell2010l) . 
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6 SUMMARY AND CONCLUSIONS 

We have proposed a new scenario for the formation of the 
warped stellar discs of young stars at the Galactic Centre 
(GC). In contrast to the work published so far, our study 
is based on the assumption that the star formation at the 
GC took place in an accretion disc but was delayed until 
the accretion disc around SgrA* was warped. To test the 
plausibility of this idea we addressed the issue of accretion 
disc warping considering two mechanisms, the Pringle (ra- 
diation) instability, and the Bardeen-Petterson effect. From 
simple arguments we showed that warping of the GC disc 
is plausible if the star formation efficiency is high, esF ^ 1, 
and the viscosity parameter a ~ 0.1. 

Our simulations are based on the idea that, some time 
after the disc became warped, the activity at the GC sub- 
sided. Rapid cooling of the disc led to fragmentation and 
star formation, while the rest of the gas was lost. Subse- 
quently, the disk was subject only to gravitational forces 
causing it to precess. In addition to the torques from the 
self- gravity of the disc, the torques induced by a surround- 
ing non-spherical star cluster were also taken into account in 
our models. We examined the long term temporal behavior 
of the warped disc models with different disc-to-black hole 
mass ratios (mass fractions), and investigated the interplay 
between the self-gravity and the cluster torques. We also 
made a comparison of our models with the observations of 
the warped GC discs. Our findings can be summarized as 
follows: 

i) When the disc mass fraction is large, Md/Mbh ~ 0.1, 
the disc breaks into distinct pieces which then precess in- 
dependently. In this case, the self-gravity torques dominate 
over the torques from the star cluster, for cluster parame- 
ters representative for the GC. The warp profiles of these 
models evolve in an oscillatory manner which is very well 
behaved for the self-gravity model. The maximum mutual 
inclinations attained (about 40°), and the time scales by 
which the disc acquires this broken configuration are com- 
parable in both models, despite a slightly faster precession 
of the cluster model. 

ii) Moderate mass discs with Md/Afbh ~ 0.01 evolve in 
two distinct ways, depending on the absence or presence of 
the torques from the star cluster. In the former case, the evo- 
lution of the model mimics that of its higher mass counter- 
part. Both the precession and the break-up of the disc occur 
on time scales which are almost exactly 10 times longer than 
for the high mass model. The model including the effect of 
the star cluster for this mass fraction exhibits a more disor- 
dered evolution. The disc breaks into parts more slowly, but 
once it arrives in a broken configuration, its pieces cannot 
stay intact due to the decreased self-gravity, and the disc as 
a whole dissolves. 

Hi) Low mass discs for which Md/Mbh ~ 0.001 precess 
on time scales of order a few 10^ yr when the disc is subject 
only to its own self-gravity, while the presence of the star 
cluster fastens the precession of the disc by about a factor of 
fourty for the estimated quadrupole moment of the star clus- 
ter. The overall inclination of the disc is not grossly altered 
during the simulations for this mass fraction. When only 
the self-gravity of the disc is considered, the disc precesses 
as a single body throughout its evolution, while inclusion of 



the torques from the star cluster cause the rings to precess 
apart. 

iv) None of our models lead to a configuration where 
parts of a broken disc would counter-rotate on the plane of 
the sky, i.e. give rise to two mutually counterrotating discs. 
A comparison of our models with the better-defined clock- 
wise rotating disc shows that the lowest mass self-gravity 
only model matches the data best. In this case in order 
for the model to mimic the observed orbital configuration 
of the stellar disc, its initial inclination determined by the 
gas disc warping mechanism has to be close to the observed 
one. On the other hand when the star cluster torques are 
taken into account, the overall shape of the warp can not be 
reproduced. The disagreement of the model including the 
star cluster might mean either that the old star cluster at 
the Galactic Centre is almost spherical so that it exerts no 
torque on the discs, or that the mass of the cluster within the 
radius of the disc is less than the value assumed here based 
on the most recent observations. Future observations of the 
old star cluster will therefore help to better understand the 
behavior of the stellar discs. 
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